Interacting electrons in disordered potentials: Conductance versus persistent currents 
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An expression for the conductance of interacting electrons in the diffusive regime as a function of 
the ensemble averaged persistent current and the compressibility of the system is presented. This 
expression involves only ground-state properties of the system. The different dependencies of the 
conductance and persistent current on the electron-electron interaction strength becomes apparent. 
The conductance and persistent current of a small system of interacting electrons are calculated 
numerically and their variation with the strength of the interaction is compared. It is found that 
while the persistent current is enhanced by interactions, the conductance is suppressed. 
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There has been much recent interest in the physics of 
interacting electrons in disordered systems p]-p^. Part 
of this attention is motivated by the large amplitudes of 
per sistent currents observed in mesoscopic metallic rings 
[po| , pT| . These values are larger by up to two orders of 
magnitude than theoretical predictions based on the sin- 
gle electron picture using the value of the mean free path 
as measured by transport experiments. Another motiva- 
tion has to do with the rich physics contained therein. 
The metal-insulator transition can be triggered by two 
different physical mechanisms: electron-electron (e-e) in- 
teractions (generally referred to as the Mott-Hubbard 
transition) and disorder (known as the Anderson tran- 
sition). Although much effort has been devoted to the 
investigation of the interplay between the two, the prob- 
lem of metal-insulator transition in the presence of disor- 
der and interactions is not yet completely settled P3,p3[ . 



Theoretically, it has been established that due to e-e 
interactions the amplitude of the persistent current (at 
zero-temperature) may be enhanced compared to its non- 
interacting value. The precise nature of this interaction 
induced modification depends on the model used and on 
the specific domains in parameter space. For spinless 
electrons in one-dimensional (ID) continuum models the 
amplitude can reach its disorder-free value for strong in- 
teractions 1^ . On the other hand for spinless electrons in 
ID lattice models a negligible enhancement of the ampli- 
tude occurs and that happens only for weak interactions 
in the localized regime pl Jl^Jl^ . When spin is taken into 
account, a sizable enhancement of the amplitude is found 
[0,0. Large enhancements occur also for 2D and 3D 
spinless electrons in lattice models for weak and medium 
ranges of interaction strengths p6 19|. 



Thus one may conclude that it is conceivable that sig- 
nificant enhancement of persistent currents may result 
in calculations for realistic 3D lattice models which take 
spin into account. Nevertheless, there still remain sev- 
eral important questions which have not yet been fully 



answered. The first, and perhaps the most interesting 
one from a general point of view, is why does e-e inter- 
action play such an important role in the determination 
of the persistent current while it is believed to play no 
important role in determining the value of the conduc- 
tance. Or, to pose the question in another way, it is 
conceivable that for the models considered, e-e interac- 
tions enhance also the conductance and therefore do not 
explain the discrepancy between theory and experiment. 
Other questions are connected with the precise origin and 
nature of the enhancement. 

In this letter we shall concentrate on the first question. 
A new expression for the calculation of the dissipative 
conductance for a system of interacting electrons is pre- 
sented. It can be written as the derivative of the per- 
sistent current at zero flux multiplied by the compress- 
ibility of the system. It is argued that in the diffusive 
regime, the derivative of the current is of the same or- 
der of magnitude its amplitude, (which is enhanced by 
e-e interaction), while the compressibility is suppressed. 
Therefore, in the same regime of disorder and interaction 
strength for which the persistent current is enhanced, the 
conductance might behave in quite a different way. The 
new formulation is then applied in the numerical evalu- 
ation of the conductance for 2D spinless electrons on a 
lattice which is known to exhibit large enhancement of 
the persistent current. The results are compared with 
the conductance as calculated via the Kubo-Greenwood 
formula (suitably adopted for interacting systems). For 
both methods of calculations, the conductance of the sys- 
tem is suppressed by the interactions, thus supporting 
the suggestion that e-e interactions might explain the 
discrepancy between theory and experiment. 

As the starting point we shall use the Akkermans- 
Montambaux p4| definition of the conductance which is 
based on the response of a system to a change in its 
boundary condition 
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where the boundary condition on the wave-function of 
the system is given by ip{xi,yi, zi;. . .■,Xj,yj, zf, . . .) = 
ip{xi,yi, zi; . . . ; Xj + L, y^, z^; . . .)e**. This form of 
boundary condition is similar to the one arising in a 
ring encompassing a magnetic flux 0. In that case 
$ = 2TT(j)/(j>Q, where 0o is the quantum flux unit. N{^, $) 
is the number of particles for a given realization of the 
disorder, with specific values of the chemical potential 
^ and Here {SN^i^i,^)) = {N^{n,<l>)) ~ (iV(^,$))2, 
and (. . .) represents an average over disorder realizations. 
Since this expression describes the conductance at a fixed 
chemical potential it is especially appropriate for a sys- 
tem coupled with the external world, for example by 
leads. It is applicable in the metallic regime in which 
the usual diagrammatic expansion is valid. In the meso- 
scopic regime, it is the case for which the level broadening 
7 is larger than the averaged single-particle level spacing 
A |p5| , which is compatible with the experimental situ- 
ation. Thus, one must be careful in applying expression 
(1) in the deep quantum limit in which the level broad- 
ening is smaller than the level separation (see Ref. pGf). 
The above definition remains valid also for interacting 
particles described by a Fermi liquid picture. A connec- 
tion between the relation (1) for the conductance and 
the Thouless formula gc = (|9^£^„/9$^|)$=o/A (where 
En is the energy of a single electron level in the vicin- 
ity of the Fermi energy) in the absence of e-e interaction 
was established analytically |2^. From numerical stud- 
ies ||2^ it seems that gd oc gc holds (for varying disorder 
strength) even in the deep quantum limit although the 
proportionality factor changes. 

We would like to express gd in terms of the persistent 
current. Altshuler, Gefen and Imry pTj have shown 
that the fluctuations in the number of particles in the 
grand canonical ensemble is connected to the disordered- 
averaged persistent current in the canonical ensemble 
(i.e., an average over different realizations of disorder 
with a fixed number of electrons A^o) in the following 
way M 
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where {dN/djj) is the averaged compressibility, and (^} is 
the averaged chemical potential for which = (/i)) = 
Nq. This connection is general and valid for interacting 
systems as well. 

Combining Eq. (If]) and (0) one obtains 
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Thus, in the metallic regime where the usual diagram- 
matic expansion is valid (7 ^ A) one can relate the 



derivative of the averaged persistent current at zero flux 
for a canonical ensemble to the dissipative conductance at 
a given chemical potential. This connection remains valid 
also for interacting systems under the previously men- 
tioned restrictions. Expression (3) is extremely useful for 
numerical calculations of the dissipative conductance at 
zero temperature for such systems since it involves only 
the ground-state properties (energy and compressibility) 
of a system with fixed number of particles. 

The definition (Eq. (|)) agrees with our physical con- 
cept of the conductance especially for systems in interac- 
tion. As pointed out by Lee , the conductivity for the 
interacting case can be written as cr = {dN/dfi)D/L'^, 
(where D is the diffusion constant and d is the dimension- 
ality of the system) which is simply the Einstein relation. 
Thus the conductance is, g = L'^^^a — {dN / dii)D / L'^ , 
which is exactly the content of Eq. This is eas- 

ily verified in the non-interacting limit where the deriva- 
tive of the averaged persistent current was calculated by 
Altshuler, Gefen and Imry j2^, and was shown to be 
0o(9/(*)/9$)jVo,$=o = (i?/L2)(A/7). For 7 = A the 
continuous spectrum conductance is recovered, and for 
7 > A the Drude formula for conductance with inelastic 
scattering is obtained. 

It is also possible to directly connect gd to the ampli- 
tude of the persistent current for the interacting case. 
In the diffusive regime, the average persistent current 
is determined by the first few harmonics of the cur- 
rent [2| j27|j30t| . The situation changes in the presence 
of interactions where it was shown analytically that the 
first harmonic describes very well the current for any 
value of the flux Q. Later on it has also been con- 
firmed numerically [^,0. Thus, for interacting elec- 
trons {I{'^))no ^ = 7r/2))Aro sin(<i>), which results 
in {dI{^)/d^)No,<i>=() ^ = 7r/2))7Vo- Inserting this 
relation into Eq. (||) one obtains gd ^ {{dN/dfi)){I{^ = 
Tr/2))]\[i^. Therefore, the conductance is proportional to 
the persistent current multiplied by the compressibility. 
This implies that strong correlations might influence per- 
sistent currents and conductance in an opposite way. 
While persistent currents are enhanced by e-e interac- 
tions, the conductance which is the persistent current 
multiplied by the compressibility (a decreasing function 
of interaction) might be suppressed at higher values of 
interactions. 

We shall now illustrate our arguments by calculating the 
conductance for a system of interacting electrons on a 
2D cylinder of circumference and height Ly. In this 
model large enhancement of the persistent current in the 
diffusive regime has been found [l^ . The model Hamil- 
tonian is given by: 

H = '^ek^jal jQk^j - Vy^(exp(z^s/£-c)Qi.,+iQfcj + h.c) 
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FIG. 1. The derivative of the averaged persistent current 
at zero flux for a fixed number of particles as a function of 
the interaction strength (in units of V). The derivatives were 
averaged over 4500 realizations. In the inset the averaged 
persistent current at a given flux for the same system is pre- 
sented. 

where j is the fermionic creation operator, ekj is the 
energy of a site {k,j), which is chosen randomly between 
—W/2 and W/2 with uniform probabihty, F is a constant 
hopping matrix element and s is the lattice constant. 
The distance \rk,j — n,p| = (min{(fc — /)^, (Lx/s — {k — 
0)^}+min{(j-p)2, [Ly/s-{j-p)f}f/'^. The interaction 
term represents a Coulomb interaction between electrons 
confined to a 2D cylinder embedded in a 3D space with 
£c = e^/s. 

Let us start by presenting numerical results for g^- Since 
we are concerned with an exact diagonalization of the 
Hamiltonian for an interacting electrons the size of the 
system is naturally limited. We consider a 4 x 4 lattice 
with m = 16 sites which in the half-filled case {Nq = 8) 
corresponds to a 12870 x 12870 matrix. As previously 
mentioned, the main interest lies in the diffusive regime. 
Therefore we must chose the disorder strength W accord- 
ingly, i.e., ^ > Lx,Ly 3> £, where ^ is the localization 
length and £ is the mean free path. We take W = 81^ for 
which ^ = 8.4s (estimated using the participation ratio) 
and £ = 0.97s (estimated from the Thouless conductance 
(7c)- We also checked that the single-electron level spacing 
distribution is close to the Gaussian ensemble prediction, 
thus confirming that the system is in the metallic regime 
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For systems with interacting electrons it is not possible 
to calculate directly the Thouless conductance, since the 
single electron energy levels are not defined. On the 
other hand, the basic ingredients needed for the calcu- 
lation of gd are easily available once the ground-state 
energy of the many-particle system as a function of flux 
is obtained by diagonalizing the Hamiltonian (^). The 
persistent current can be calculated via the well known 
relation (/($))jVo = -(27r/0o)(a^;($)gs,JVo/a$), (where 
E{^)gs,Na is the ground state energy of an interacting 
system of A^o particles), and (dfi/dN) = (Egs^Ng+i ^ 
2Egs^No + Egs^No-i)- Thus, once the ground states en- 
ergies for systems with Nq ± 1 particles are known, the 
conductance gd can be immediately calculated. 



FIG. 2. The dissipative conductance gd as a function of the 
interaction strength. In the inset the inverse of the averaged 
compressibility as function of the interaction strength for the 
same system is presented. 

0.40 I ■ 1 ■ 1 ^, 1 : 1 




O.OO I ' ' ' ' ' ' ' 1 

O.OO o.o I 0.02 0.03 0.04 

FIG. 3. The correspondence between the Kubo conduc- 
tance Qk an the dissipative conductance gd- For strong inter- 
actions (ec > 'iVjgd < 0.015) a good fit with a linear relation 
gk ~ 13.9gd + 0.04 is obtained. In the inset, gk as a function 
of the interaction strength is presented. 

The derivative of the persistent current at zero flux for 
Nq — 8 averaged over 4500 realizations is presented in 
Fig (|^). It can be seen that an enhancement of the 
derivative as function of the interaction strength Sc is 
obtained. This is similar to the enhancement of the cur- 
rent at $ = 7r/2 shown in the inset. Thus our assumption 
((9/($)/9$)Aro,<[.=o ~ = 7r/2))wo is validated. 

The inverse of the compressibility at Aq = 8 averaged 
over 500 samples is shown in the inset of Fig (^ . It can 
be seen that for Sc > 2V (dfi/dN) ~ st^^. The dissipa- 
tive conductance is presented in Fig (|^). For small values 
of Ec, gd shows a substantial decrease since the deriva- 
tive of the current is only weakly enhanced and the main 
influence on the conductance comes from the compress- 
ibility. For larger values of interactions {Sc > 2V) the 
enhancement of the current is compensated by the com- 
pressibility and the conductance slowly decreases. Thus, 
for all values of Ec the conductance is not enhanced by 
the interactions. This strongly points toward e-e interac- 
tions as a possible explanation for the large amplitudes 
of persistent current measured in experiments |^,^ . 

A useful check of our formalism is to compare the val- 
ues of gd to the values obtained from a many-particle 
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formulation of the Kubo-Greenwood conductance gk for 
the same system. Following Kohn [ p3|j3^ ] the real 
part of the conductance may be written as: gk = 

(8^ViV)EJ{«l-/.|0)pe„,o7«o + 7')-', where |0) 
is the many-particle ground-state, Jx is the current op- 
erator, and £0,^0 = £a — Eo- This is a very cumbersome 
calculation since it involves calculating the many-particle 
low-lying eigenvalues and eigenvectors for each realiza- 
tion of disorder. We chose the inelastic broadening to 
be of the same order as the single electron level separa- 
tion, i.e., 7 = 0.7V. The results for gk averaged over 
180 realizations are plotted in the inset of Fig. ^ As 
in the relationship between g^ and g^ [ p4|j2^ g^. is an 
order of magnitude larger than gd. gk seems to follow 
gd for Sc > '2,V (for which the excitation separation is 
much bigger than 7), while there is no sharp decrease 
for small values of interactions. This may be clearly seen 
in Fig. m where gk is plotted as function of gd. For 
Ec > 2V a clear linear relation gk = 13. 9^^ + 0.04 is 
obtained, which can be compared to the non-interacting 
relation gk = 8.9^^ + 0.04. For Ec < 2V gk is almost 
constant while gd decreases. This might be connected to 
the transition of the statistical properties of the many- 
particle energy levels for weak interactions |35| and to 
the influence of interactions on the inelastic broadening 
7 currently under investigation. Nevertheless, in both 
methods of calculation the conductance never increases 
as function of interaction strength. 

In conclusion, we have presented a new method for calcu- 
lating the conductance of an interacting electronic system 
in the presence of static disorder. This formulation clari- 
fies the difference in the interaction dependence between 
the persistent current and the conductance. Since the 
compressibility decreases as function of e-e interaction, 
the conductance is always less enhanced (or probably 
suppressed) compared with the persistent current. Thus, 
e-e interactions are a probable candidate for explaining 
the discrepancy between theory and experiment. This 
method has also the advantage of being dependent only 
on ground state energies of neighboring systems, which 
are much easier to calculate numerically than the full 
spectrum of eigenvalues and eigenvectors. 
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